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Chapter 1 
Introduction 


The electromagnetic behavior of an ordinary isotropic, linear, homogeneous material 
is completely specified by the complex constitutive parameters, or the permittivity 
and the permeability, of the material. The permittivity and the permeability of 
the material are denoted as t and [i. If these constitutive material parameters are 
unknown, they must be determined before any electromagnetic problem involving 
the material can be solved. 

The literature is rich with a variety of techniques to determine these constitutive 
parameters [1,2]. Of these techniques, one appears to be very practical to determine 
the permittivity of a non-magnetic material which would naturally be in a slab form. 
This technique relates the reflection coefficient of an aperture in a material coated 
groundplane. This technique is not new and was first used by [4, 5] and by any other 
in recent times [6, 7, 8]. 

Unlike the closed form parameter determination using reflection and transmis- 
sion waveguide measurements [9], a numerical search has to be performed between 
measured and calculated values for this aperture technique. A Newton- Raphson it- 
eration scheme is often used using the complex reflection coefficient if the material 
is lossy. For cases when the material is lossless, the magnitude or the phase of the 
reflection coefficient is sufficient to make a parameter determination. 
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Material slab 


The parameter determination using an aperture technique as several advantages 
over other techniques. This technique is convenient for planar samples since minimal 
sample preparation is required and is non-destructive. Errors resulting from poor 
sample shape (waveguide air gaps) are eliminated Additionally) the phase reference 
plane is exactly determined by the groundplane surface unlike in waveguide mea- 
surements where the possibility of some shift always exists. Phase reference stability 
is very important in the determination of large permittivity values. 

The work in this report expands the technique to include additional reflection 
coefficient expressions for different material configurations and apertures. The aper- 
tures of intersect are the rectangular and coaxial apertures as shown in Figures 1.1 
and 1.2. The coaxial aperture is of interest due to the potential of broader bandwidth 
measurements, especially at lower frequencies. It should be noted that a commercial 
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Ground plane 



Material slab 


Figure 1.2: Coaxial aperture 


coaxial probe is available but the software often used with it does not account for 
sample thickness [10]. 

The material configurations presented and studied in this report are shown in 
Figures 1.3 and 1.4. This material configurations allows for the study of material- 
groundplane air gaps, PEC (perfect electric conductor) backed samples, sheet re- 
sistivity determination and the equivalence of sheet resistivity modeling with lossy 
material layers of finite thickness. 


3 





Resistive 

Sheet 



Resistive 

Sheet 


Figure 1.3: Rectangular aperture 



Resistive 

Sheet 



Figure 1.4: Coaxial aperture 
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Chapter 2 

Numerical Model 


The numerical model upon which the permittivity de-embedding procedure is based 
is developed in this chapter. The functionals in the de-embedding procedure are 
based on the reflection coefficient of the rectangular or coaxial aperture, covered by 
a dielectric slab/resistive or PEC sheet configuration. The reflection coefficient is 
computed from the aperture admittance. The derivation of the aperture admittance 
for both aperture shapes is based upon the rectangular aperture work of Rudduck 
[12]. In each of the two sections to follow, the basic derivation of the rectangular 
or coaxial aperture admittance is presented, followed by the specifics for the various 
dielectric slab/resistive or PEC sheet configurations under investigation. 

2.1 Rectangular Aperture 

The admittance of a rectangular aperture is developed here. The analysis involves 
representing the fields external to the aperture in terms of the Fourier integral of 
their plane wave spectra. The aperture admittance is obtained by using a variational 
assumption for the aperture total field distribution and matching the boundary 
conditions for the various dielectric slab/resistive or PEC sheet configurations under 
study. 
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The analysis starts with the representation of the electric vector potential F in 
terms of the rectangular electric scalar potentials 0 and ip: 

Fo, 1,2 = 00,1,2* + 00,1 l2 y. ( 2 * 1 ) 

Referring to Figure 2.1, it can be seen that the potentials are functions of x, y , 
and z, defined only for z > 0. The subscript 1 refers to Region 1, the dielectric 
slab in contact with the groundplane. The subscript 2 refers to Region 2, a second 
dielectric slab which may be on the first dielectric slab. The free space region above 
the dielectric slab(s) is denoted as Region 0 and all quantities associated with that 
region possess the subscript 0. The scalar potentials are solutions to the scalar wave 
equation 


(v J + *s, u ) ( } , 

( 00,1,2 J 

(2.2) 

where the wavenumbers fco,i, 2 are fP ven by 


k>l = u> 2 n 0 e 0 , 

(2.3) 

for Region 0, and 


k 2 2 = u; 2 /io£i,2 = ^o Cr i,2 1 

(2.4) 

for Regions 1 and 2. In general the relative permittivities e r , , 

are complex, repre- 


senting possible losses in the slabs. 

The scalar potentials can be expressed in terms of their plane-wave spectra. 
In Regions 1 and 2, the total spectra are comprised of the incident and reflected 
spectra, respectively denoted by I and R. In Region 0, the total spectra are each 
comprised of one transmitted spectrum, denoted by T. Thus the plane wave spectral 
representation of the scalar potentials in Region 1 are given by 

*(.,»,*) = [/*,(*. AK'*"' + A)* +, ‘"‘] 

e~* klX e~* kyV d.k x dk v (2.5) 
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GROUND PLANE 




REGION 1 (e lf /4 0 ) 

POSSIBLE RESISTIVE 
OR PEC SHEET 

REGION 0 (e 0< Aio) 



Figure 2.1: Dielectric slab/resistive sheet configurations on rectangular aperture 
groundplane. 
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and 


e -i k ** e -i k »«dk x dk y ] 


( 2 . 6 ) 


in Region 2 by 

V’at*,!/,*) 


irj [/*(*., + «*(*•• 


(2tt) 

vdk x dk v 


(2.7) 


and 


</»2 (*,!/,*) = 


(Sj5 jC II M*.. + «*>(*»> t ») e+ ’*‘’*] 

e~^ ktX e~ 3kyV dk x dk y \ 


( 2 . 8 ) 


and in Region 0 by 

*,(»,»,*) = T *>(* (2-9) 


and 


*>(*,»,«) = Tjk /“ /” T i ,(k„K)e-‘ k -‘-e-> l -’e-‘ k "‘ik,ik,. (2.10) 


(*»)■ 

The propagation constants are given by 


fc*0,l,2 — i\Z^O,l,2 — 


( 2 . 11 ) 


where the sign is chosen so that the radiation condition is satisfied. Thus the sign 
for each propagation constant is chosen so that 


»(*,,,)>« 


( 2 . 12 ) 


and 


*(** 0 ,!,,) < 0 . 


(2.13) 
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From vector potential theory, the electric and magnetic fields are related to the 
electric vector potential by 

E = -V x F ( 2 . 14 ) 

and 

H = K'-’ F + V(V ' F >1 ■ (2-15) 

Therefore the rectangular components of the electric and magnetic fields in Region 


1 are given by 



E X1 (® > j/i 2 ) = 

1 [OO AOC p 

(2tt) 2 I-cc I - 00 6 1 + i kzi ^ ** “ ] 



e~*** e ~ ikyV dk x dk v , 

(2.16) 

E vA x ,y,z) = 

\ /*oo *00 _ 

(2tt) 2 Loo Loo^ kziI ^ e J 1 - 11 ] 



e-^ x e-^«dk x dk y , 

(2.17) 

H Xl (x,y,z) = 

1 roc foe r r.2 

(i-xv! I 1 — £ K e ~ 3k ‘ i z + e3kH 1 - 

{2ry j -00 J -00 { jujfio L ™ J 

— [/* e-^i * + *1 } e-’ k **e-J k v»dk x dk v , 

JUJflQ J J 

(2.18) 

and 



H vi ( x ,y,z) = 

jC/o \ h ' C ~“"' * + ** e'*- 1 *]} e-*“ iKik,-, 

(2.19) 

in Region 2 by 



E*,( x >y, z ) = 

(2t) 2 I-oc Loo 3 2 *’ ] 



e- jk * x e- jk * v dk x dk tn 

(2.20) 


(2tt) 2 L 00 I -00 t 7 **’ 7 ^ 6 J 2 *2 ] 



e~ 3kxZ e~ jkyV dk x dk y , 

(2.21) 
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e~^ x e-^dk x dky, 


ju>Ho 


and 


e- ik * x e- ik ">dk x dky-, 


k X fcy 
jw/io 


and in Region 0 by 


*"(*■»■*)= (SF /I /I 


E, 


VO 


(*,»,*) = 7 TvT T iKT^-i^e-i’-e-^ik.ik,, 

yZ 7T J J —oc J — oo 


e -i k *o z e ~i k * x e ^ v dk x dk y , 


and 


jfftf °( a; ’ y ’ 2 ) (2ir) 2 J-oJ- 


00 f “ r An - . 


krky 


0 "y J _ 

(27r) 2 J — 00 y-00 [ jwfto jupo 

e~ jk '° * e - jk * x e- jk '> v dk x dk v . 


4>o 


( 2 . 22 ) 


(2.23) 


(2.24) 


(2.25) 


(2.26) 


(2.27) 


Note that in the above the fc*, fc„ functional dependence is dropped from the spectra 
and are understood. 

The solution to the fields can be obtained by satisfying the boundary condi- 
tions at the groundplane/ dielectric slab interface and other interfaces present in the 
dielectric slab /resistive sheet configurations if the total fields in the aperture are 
known. However, the total fields in the aperture are not known, since they consist 
of an infinite sum of waveguide modes with unknown coefficients. In order to find 
the coefficients, a Moment Method solution can be set up in a manner similar to 
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that which has been performed for the coaxial aperture case [11]. For the purposes 
of this investigation a simple variational solution will be used [3]. The variational 
solution consists of the assumption of a form for the aperture electric field whereby 
first-order deviations of the assumed field from the exact field will be manifested 
as second-order errors in the admittance. The electric field is assumed to be of the 
form of the dominant TEio mode. This aperture field assumption works well in 
the case where the waveguide dimensions are such that all higher-order modes are 
evanescent [12]. In such a case any deviations of calculated admittance values from 
measurements are attributed to the exclusion of higher-order modes. 

Thus the total electric field in the plane z = 0 is assumed to be 


E(*»y»0) = e 0 = < 


& cos f* 


0 


in the aperture 
outside the aperture 


(2.28) 


where the factor normalizes the electric field so that 

/ OC i»00 

/ E(x,y,0) • E(x,y,0)dxdy = 1. 

-OO «/ — OC 


(2.29) 


Referring back to Figure 2.1, it can be seen that b is the length of the longer edge 
of the waveguide, and a is the length of the shorter waveguide edge. 

The rectangular waveguide aperture admittance, Y a , is given by 
Ya = HJSo jgfeM x H(g,y, 0)dxdy 


^ oo roc 

/ / (®j 2 /j 0) Hy x (x, j/j 0}dxdy. 

J — oc J — OO 


(2.30) 


Evaluation of (2.30) is accomplished through the use of Parseval’s Theorem in two 
dimensions in a rectangular coordinate system: 

/-oo I- oc E *' H «' dxd y = J- oo I- oc e *i‘ H » dk * dk *> ( 2 * 31 ) 

where £ Xl and 7i Vi are the Fourier transforms of E Xl and H Vl , at z = 0. These 
transforms can be found from (2.16) and (2.19) and are given by 


£ri (^IJ fcj/) — j^Zl + ■ftt/’l] 


(2.32) 
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and 


W yi {k x ,k v ) — . [(fcj & y )(Ah d" ** ) + -^>1)] • (2.33) 

JW/io L 

Thus the aperture admittance can be written as 

y ° = ttvT r 

yZ'K J J — oo J — oo 

{i [<*> - W* +*>.)- W* + «*.)]} < 2 - 34 ) 

The unknown quantities in (2.34) are the incident and reflected spectra of Region 
1. The boundary condition of continuous tangential electric field at z — 0 for each 
pair of spectral values (k x ,k y ) gives the relationships between the Region 1 incident 
and reflected spectra as 

/2fcsin (^) cos (^) 

-jK [-4i ~R*i] = 4ir \/ 7 kjy _ *252) = f ( 2 - 35 ) 


j^Zl [I<fri R<t>i ~ 


(2.36) 


where the Fourier transform of the aperture electric field is denoted by / for the 
sake of brevity. The substitution of (2.35) and (2.36) into (2.34) yields 

= (SJi £ £ [M- 1 ^ ' W] } f**.. ( 2 . 37 ) 

In order to make the numerical evaluation of the admittance more convenient, the 


following change of variables is made: 
k x = fco/3 cos a 


(2.38) 


k v = Aso/3 sin a. 

Therefore the aperture admittance can be written as 


—2 fi 2 sin a cos ak^I^ j ifkafidadfi, 


(2.39) 


(2.40) 
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where Y 0 is the admittance of free space. In order to use the aperture a dmi ttance in 
the calculation of the reflection coefficient, it must first be normalized with respect 
to the TE 10 mode waveguide admittance: 


Y 

y r — _ — 

n — — 


Y. 


(2.41) 


where \ g is the wavelength in the waveguide. Thus, the reflection coefficient of the 
rectangular waveguide aperture is given by 
1-K 


r = 


1 + Y n 


(2.42) 


A few items about the integration in (2.40) must be stated. The integration 
is performed along the path in the complex plane shown in Figure 2.2. The 
Imjj}} 



Figure 2.2: Integration path in the complex /3 plane, with surface wave poles and 
parallel plate waveguide pole (PEC sheet backed cases only). 


integration path is deformed from the real axis so as to avoid surface wave poles. 
Surface wave poles can occur on the real axis when there is at least one lossless 
dielectric layer in any of the five configurations of Figure 2.1. These poles occur in the 
interval 1 < (3 <e r , where e r is the greater of the two real relative permittivities, or 
is the only real relative permittivity. Loss in the dielectric layer(s) causes the layer(s) 
surface wave poles to move below the real axis. The integration path starts at 0 = 
0, goes to 1.0 -f- j'0.4, then to |e rg | -|- jOA, 2|e rff |, and then oo -f jO. This integration 
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path was determined experimentally and found to work well. The integration path 
crosses no poles during its deformation from the real axis, so no residues need to be 
computed and included in the integral. 

If there is a PEC sheet present, an additional complication occurs when the 
frequency is high enough to excite one or more parallel plate waveguide modes. A 
mode manifests itself as a pole in the complex plane by appearing at /3 = 0 at 
the cutoff frequency, and then moves along the positive real axis as the frequency 
increases. No provision has been included in the codes for this contingency, since 
the codes are to be used only at frequencies where there are no propagating parallel 
plate waveguide modes. 

However, the solution for the reflection coefficient is by no means complete. The 
aperture admittance is still dependent on the unknown incident spectra and I ^ , 
which are dependent on the boundary conditions given by the particular dielectric 
slab/resistive or PEC sheet configuration in contact with the aperture ground plane. 
The following subsections give the derivation of 1^ and I ( j >l for the five configurations 
under investigation. 

2.1.1 Two Dielectric Layers - Free Space Backed 

Figure 2.3 shows two dielectric layers, free space backed, on a rectangular aperture 
groundplane. The width of the first layer is d u and the width of the second layer 
is d 2 - d\. The incident spectra and are determined by the enforcement of 
the boundary conditions at the three interfaces at z = 0, d\, and d 2 . The bound- 
ary condition enforcement occurs at each set (k x ,ky) of spectral frequencies. The 
enforcement of boundary conditions results in a system of linear equations with the 
spectra as unknowns to be solved. The boundary conditions with their resulting 
equations are listed here: 
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Figure 2.3: Two dielectric layers, free space backed on a rectangular aperture ground- 
plane. 

The boundary condition of continuous tangential electric field across the interface 
at z = 0 gives the following two equations: 

-jKh, +jk Zi R^ = /, (2.43) 

J k z ihi ~ jK R^n = 0. (2.44) 

The boundary condition of continuous tangential electric field across the interface 
at z = d\ gives the following two equations: 

-iKh^ ikzidl + jk, l R* l e’ k 'i dl = -jk z J^e-’ k ‘* d ' -f jk Z7 e* kx * dl , (2.45) 
j k n hi e~* kt i dl — jk Xl R i p 1 e 7 **! dl = jk Z2 It 2 e~ ikz ' dl - jk Z3 R <h e ik '> d ' . (2.46) 

The boundary condition of continuous tangential magnetic field across the interface 
at z = di gives the following two equations: 

(*? - k l) [hi e~ 3k ‘ l dl + R*i e*t ] - k x k v [j 0l e~* kz ' d ' + e*i d ' ] = 

(*2 ~ k l) [l^ kz ^ + - k x k y \l^e~^ d ' + R^e’^] , (2.47) 
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(J fc? - fcj) + R+ l e ik '' dl ] = 

(k 2 2 - k 2 y ) [l^e-’ k ‘> dl + i^e i ^ <<1 ] - k x ky [l <h e- jk ‘* dl + R^e }k ^] . ( 2 . 48 ) 

The boundary condition of continuous tangential electric field across the interface 
at z — di gives the following two equations: 

-jk z J^e- jk '* d ' +jk„R^e jk '’ d * = -jKT^e~ ik ^ d \ ( 2 . 49 ) 

jk^e-^-jk^^ = jKT^e-*^. ( 2 . 50 ) 

The boundary condition of continuous tangential magnetic field across the interface 
at z — d 2 gives the following two equations: 

{k 2 2 - k 2 x ) [l^e-’ k " d ' + R^e 3 ^] - k x k v [l^e~ jk ^ + R^e jk ^\ = 

{k 2 0 - k 2 x )T <t>0 e- jk '° d > - KKT+e-*** , ( 2 . 51 ) 

(fc 2 2 - k 2 y ) [l^e~ 3k ^ + R^e***] - Kk v [l <h e~ ikz ^ + R^e***] = 

(kg - kl)T* 0 e-*'° d ' - k x k y T^e- 3k '° d > . ( 2 . 52 ) 

These ten linear equations in the ten unknown spectra form this matrix equation: 
Z\\ 0 — Z\\ 0 0 0 0 0 0 0 7^ ’l 

0 — Z\\ 0 Zn o 0 0 0 0 0 7 ^ 

z 3i 0 Z33 0 Z35 0 Z37 0 0 0 

0 z 3 i 0 z 3 3 0 z 35 0 z 37 0 0 R4 , 

z 51 z 52 Z53 Z54 z 55 z 56 Z57 z 58 0 0 7^2 

Z 61 Z51 Z 63 -253 2 6 5 Z55 Z 6 7 2 57 0 0 I<fa 

0 0 0 0 Z75 0 Z77 0 Z79 0 R^ 

0 0 0 0 0 Z75 0 Z77 0 Z79 

0 0 0 0 z 95 Z96 Z97 Z98 Z99 29^ T^ 0 

0 0 0 0 Z_45 Z95 Zj 17 Z97 ZA9 299 Tfo 
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where 


Zll 

— 

~jk Zl , 

265 

= 

~(H \-K)e-^, 


-^31 


—jk Zl e~i k *i dl , 

267 

= 

~{k\ — k y )e + i k *2 dl , 


■2-33 

= 

jk Zl e +]k *i dl , 

275 

= 

-jk zj e- jk *2 d \ 


^35 

= 

jk^e~ jk ‘2 dl , 

277 

= 

jk Z2 e +jk '2 


-237 

= 

-jk Z2 e +jk ‘? dl , 

279 

= 

jk z 0 e- jk ‘° d *, 


-251 

= 

-k x k v e~ jk ^ d *, 

295 

= 

-k x k y e~ ik ‘* d \ 


-252 

= 

(k 2 — kl)e~* kz i dl , 

296 

— 

(kj - kl)e- jk ‘i J \ 


^53 

= 

-k x k y e +jkl i dl , 

297 

= 

—k x k v e +jkx 2 d2 , 

(2.54) 

2 54 

= 

(k 2 - kl)e+*‘ i*. 

298 


(k 2 - kl)e^ d2 , 


^55 

= 

k x k y e~i k ‘* dl , 

299 

= 

k x k y e~ ikt o d 2 , 


-256 

= 


29/1 

= 

-(*o 2 -* 2 )e-•*‘ 0 ^ 


-257 

=r 

k x k y e +]k ‘i dl , 

2/15 

= 

{k\ ~ k 2 v )e-’ k ‘* d2 , 


-258 

= 

-(fe 2 2 -Jb 2 )e + ^ <<1 , 

2/7 

= 

( k\ - k 2 y )e+*” d2 , 


2-61 

= 

(*? - k y )e~* k *i dl , 

2/9 

= 

-(k 2 0 -k 2 u )e-**o d >. 


-263 

~ 

(Jfe 1 2 -* 2 )e + '**i dl , 





Note that 

the hexadecimal number A is 

used 

in the place of the s 

ubscript 10 for the 


sake of simplicity. The matrix equation of (2.53) can be reduced into a four equation 
system in the four unknowns 7^ 2 , 1 ^ , , and R ^ . These four unknowns have a 

simple relationship to the desired unknowns 7^, and 7^ . Thus the solution for the 
spectra 7^, and 7^,, is given by 


#11 #12 #13 #14 




' V x ' 

#21 #H #23 #13 


I<fo 



-^31 #32 #33 #34 




0 

#41 #31 #43 #53 


Rfa 


0 


where 


Z\\ = k x k. 


x^y 


1 -JT^-cot (k Xl di) 

k Zl 


9 j^x 2 dl 


(2.55) 


(2.56) 
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Z\2 

— 

(*’ - *!) - ir-M - *?) “*(*., <*. )1 e-*-** . 

(2.57) 

Z\3 

= 

fcife u 1 + cot(/s z ,dx) e +3k ‘* dl , 

(2.58) 

Z\4 

= 

(K - kl)+j%Hkl - *?)cot(M.)l e*' 4 -*, 

**l 

(2.59) 

Z 21 

= 

(kl - k\) - j^-(kl - *?)«*(*„ <f,)l «-*»*, 

(2.60) 

Z 23 

= 

(kl - kb+j^ikl - k’)cot(k„ d «+*■>*, 
K *\ 

(2.61) 

Z 31 

= 

k x k y fe- - ll e~ ik " d \ 

lAo 

(2.62) 

Z 32 

= 

(kl-k^-(kl-k 2 2 )]e- 3k ^, 

k z 0 

(2.63) 

Z 33 

— 

-k x k v fe- + ll e +jk ” d \ 
.*0 

(2.64) 

Z 34 

= 

- (*J - *2)^ + (*. - *1)1 

(2.65) 

Z 41 

= 

(ki - ki)^- - - d «-*■**. 

(2.66) 

Z 43 

= 

- (K - ‘Sir 2 - + (*J - *i)l e**-’"’’ 

*z 0 

(2.67) 

Vx 

= 

k f k ‘(tot(k„i l )-j)e*’ t “ J 'f, 

“*1 

(2.68) 

v 2 

= 

^ V l (cot(/s Zl di) j)e +J * l l /, 

**1 

(2.69) 


h, = f*” 7 *"] (e» k -’ d ‘R*, - e-’'-’ d 'I*,) + |c ° t( ^ l)+ ^/ , (2.70) 

sm(* Zl ai) v ' ■‘*ii 

*■ = < 2 ' 71 > 

Note that the substitution of variables in (2.38) and (2.39) is needed in (2.56)-(2.71) 
for the purposes of computation of the aperture admittance. 
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2.1.2 Two Dielectric Layers - PEC Sheet Backed 


Figure 2.4 shows two dielectric layers, PEC sheet backed, on a rectangular aperture 
groundplane. The width of the first layer is d \ , and the width of the second layer is 
di d\. The PEC sheet is located at z = d 2 . The incident spectra I^ l and I ^ are 
determined by the enforcement of the boundary conditions at the three interfaces at 
z = 0, rfi, and d 2 . The boundary condition enforcement occurs at each set (k x , k v ) of 
spectral frequencies. The enforcement of boundary conditions results in a system of 
linear equations with the spectra as unknowns to be solved. The boundary conditions 
and resulting equations at z = 0 and z = d\ are the same as for the free space backed 
case of 2.1.1. 



Figure 2.4: Two dielectric layers, PEC sheet backed on a rectangular aperture 
groundplane. 

The boundary condition at z = d 2 is that of vanishing tangential electric field, 
which produces the two equations: 

-j k X7 1^ e~ jk ‘ 2 d * + j e jk ‘* =0, 

- jk^e^ = 
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0 . 


(2.72) 

(2.73) 


(2.43)-(2.48), (2.72), and (2.73) form a linear system of eight equations in eight 
unknown spectra, given by the matrix equation: 

zw 0 — z n 0 0 0 0 0 Ifa 

0 —zn 0 Zn 0 0 0 0 Ifa 

z 3 i 0 z 33 0 Z 35 0 Z 37 0 72^! 

0 Z31 0 Z33 0 Z35 0 Z37 R *, 

Z 51 z 52 Z53 Z54 Z55 Z56 Z57 Z58 Ifa 

Z 61 *51 ^63 2 53 *65 Z S 5 z 67 z 57 Ifa 

0 0 0 0 Z75 0 Z77 0 Ufa 


where 


L ° 


0 000 

275 

0 

Z 77 J [ J 

*11 

— 

-iK 1 , 

256 

rr: 

~i k 2 ~ k Z)e- jk ‘2 dl , 

*31 

= 

-jk Xl e'’ k 'i d ' , 

257 

= 

k x k y e +3k *i dl , 

*33 

= 

jk Zx e +jk ‘i dl , 

258 

= 

— (fcf — A£)e +, *** dl , 

*35 

= 

jk tJ e- }kl i di , 

261 

— 

(*? - ky)e~* kz i dl , 

^37 

= 

— jk z ,e +ik ** dl , 

263 

= 


*51 

= 

-k x k v e~^ d \ 

265 

= 

-(*2 - fcy)e - - ,tijdl , 

*52 

= 

(k*-k 2 x )e-*'i d \ 

267 

= 

-(kl - kl)e +3kz ^ , 

*53 

= 

-k x k v e +jk *i dl , 

275 

= 

-jk ZJ e~ jk ‘* d7 . 

*54 

— 

(% - kl)e +*«* S 

277 

= 

jk Zl e + * k ’i d7 . 

*55 

= 

k x k y e~ 3k ^ dl , 





/ 

0 

0 

0 

0 

0 

0 

0 


(2.74) 


(2.75) 


The matrix equation of (2.74) can be reduced into a two equation system in the 
two unknowns /^ 2 and Ifa. These two unknowns have a simple relationship to the 
desired unknowns Ifa and Ifa . Thus the solution for the spectra 7^, and Ifa is given 
by 

v.Zu - V 2 z 12 


I fa — 

I fa 


Zn ~ Z\ 2 Z 2 \ 
VjZrn - Vi Zn 
Zn ~ Z\ 2 Z 2 \ 


(2.76) 

(2.77) 
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where 


Z\\ = k x k y ^1 — cot^^dx)^ e jkx i dl 

+ (l +Jj a cot(*„<i 1 )j e+>»n* e -«^.* , 

Z.2 = f(*2 - *1) - j(*2 - Jtf)^cot(*„rf,)l 

+ m - kl) + j(kl - jf)Ja.cot(t.*)l 

L **i 

«»i = f(*; - *J) - J(*; - «.*(*„<*, )1 

+ (K ~ * 2 ) + i(^y — k l)j~ cot (k^di) e +J kx 2 d l e ~ 2 J kx 2 d 2 } 
L k z, 

Vi = ~T~( coi (kz 1 d 1 ) — j)e +]kt ' dl f, 

™2i 

V 2 _ (_ y — (C0t(fc n dl) ~ j)e +Jk *l dl f, 

**l 


(2.78) 

(2.79) 

(2.80) 
(2.81) 
(2.82) 


Ah 

A> 1 


■ iKJK) / 2 jk M ,d2 +jk z ,d 1 -jk„dA j [eot^diJ + j]/ 

2j sin(A: 2] di) ' + 2 fc 2l ’ (2 ' 83) 


_ (k z Jk Zi ) / 2 jk d7 + J k^d l _ -jk^ \ j 

2j m(k ai di)y ) 1<h ' 


Note that the substitution of variables in (2.38) and (2.39) is needed in (2.78)-(2.84) 
for the purposes of computation of the aperture admittance. 


2.1.3 One Dielectric Layer - Resistive Sheet Backed 

Figure 2.5 shows one dielectric layer, resistive sheet backed, on a rectangular aperture 
groundplane. The width of the single layer is d x . The resistive sheet is located at 
z — d\. The incident spectra 7^, and I ( j >l are determined by the enforcement of the 
boundary conditions at the two interfaces at z = 0 and d x . The boundary condition 
enforcement occurs at each set {k x , ky') of spectral frequencies. The enforcement 
of boundary conditions results in a system of linear equations with the spectra as 
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unknowns to be solved. The boundary conditions and resulting equations at z — 0 
and are the same as for the free space backed case of 2.1.1. 



Figure 2.5: One dielectric layer, resistive sheet backed on a rectangular aperture 
groundplane. 


The resistive boundary conditions that the electric and magnetic fields obey at 
z — d\ are 


z x (Eq — Ei) — 0 


(2.85) 


z x z x E = — Ri x (H 0 — Hi) 

(2.85) gives rise to these two equations: 

+ jfc ZlJ Rv,y fc * ldl = —jk X0 T^ 0 e ~* k *° dl , 
jk tl I^e~ ik ^ - jk ti R+ l e ik 'i d ' = 

(2.86) gives rise to these two equations: 

+ **■***] 


( 2 . 86 ) 


(2.87) 

( 2 . 88 ) 
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jufl 0 

r* 0 2 -* 2 T “x^v T 

** V*o J. . .. -‘-fai 


krk 


jk Z 0 T^e ik *° dl 


= *{.’ 


J^o 
*0 - * 2 


krk. 


JUfi o 
*? - * 2 


r 


D jkx 0 d\ 


(2.89) 


k X fay 


'x^rp "x^y m 

0 • -* *l>o 

JVfl 0 

'*e _ifc «» rf ‘ + R (t , l e }k ‘i d '] 


[< 


J^/io 


[^i e ifcil dl + 72^ e ifc *> dl ] | . (2.90) 

(2.43), (2.44), and (2.87)-(2.90) form a linear system of six equations in six unknown 
spectra, given by the matrix equation: 


where 


Zll 

0 

-Z11 

0 

0 

0 


Ah 



0 

— *n 

0 

2ll 

0 

0 


hi 


0 

*31 

0 

^33 

0 

235 

0 




0 

0 

■231 

0 

233 

0 

235 


i^i 


0 

^51 

252 

^53 

254 

255 

256 


Ah 


0 

^52 

262 

“254 

264 

256 

266 . 


** 


0 

2ll = 




254 

~ 

R 

** fc v c +i^z 1 di 

j^o 


^31 ~ 

-jk Zl e 

-J^zj 

d 1 

> 

255 

= 3 

[‘ 

+ e - 
w/io 


(2.91) 


o -jk* 0 di 


o+jfc*! d\ 


^56 

^62 


^52 — 

^53 = 


JW^io 5 

_ _ flMzMle-;**! di 

JW/i 0 ’ 

’ 64 11 jw/i 0 e » 

= -i [*„ + 

5 


■233 = jk Zl e 

■235 = jk Zo e~ }kl o d ' 

■251 = 


(2.92) 


The matrix equation of (2.91) can be reduced into a two equation system in the two 
desired unknowns 1^ and 7^, . Thus the solution for the spectra 1^ and 7^, is given 

by 

Vj Zi 2 — V2Z12 


Ah — 


Z\\Z 22 — Z \ 2 


(2.93) 
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(2.94) 


hi = 


V2Z1I - VjZl 2 
^11^22 — ^12 


where 


Z\\ = 2JZ 


( fc? - Ml 

jUflQ 


cos(k Zl d 


1) + 2 (1 + A 


w - *;) 

k Zo ufi 0 


Z u = (co8(fc Xl <ii) + 

jw^to \ **o / 


Z 22 = — 2J2^ — — cos 




<*„*)- 2 (i+« 


w - 

k Zo u>no 


Vi = 


T z Rk x k y 
v 2 — 

W/lo 


1 + ^rX~ + *, JJ 7 ’ 


^ k Zl 6in(fc Zl <fi), (2.95) 

. (2.96) 

)k l8 in(Mi)> (2-97) 

(2.98) 

(2.99) 


Note that the substitution of variables in (2.38) and (2.39) is needed in (2.95)-(2.99) 
for the purposes of computation of the aperture admittance. 


2.1.4 One Dielectric Layer - Free Space Backed 

Figure 2.6 shows one dielectric layer, free space backed, on a rectangular aperture 
groundplane. The width of the single layer is d\. The incident spectra and 1 ^ 
are determined by the enforcement of the boundary conditions at the two interfaces 
at z = 0 and d\. The boundary condition enforcement occurs at each set (k x ,k y ) of 
spectral frequencies. The enforcement of boundary conditions results in a system of 
linear equations with the spectra as unknowns to be solved. The boundary conditions 
and resulting equations at z = 0 and are the same as for the free space backed case 
of 2.1.1. The boundary conditions on the electric field at z = d\ are the same as for 
the resistive sheet backed case of 2.1.3, that of continuous tangential electric field. 

The boundary conditions on the magnetic field at z = d\ is that of continuous 
tangential magnetic field, giving the following two equations: 

(*1 2 - *x) [/*e-*“* + - k x k y [j*e-*i* + **e*-»*] = 

(*o - K)T^e-’ k ‘° d ' - k x kyT^e-^ ° dl , (2.100) 
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GROUND PLANE 



Figure 2.6: One dielectric layer, free space backed on a rectangular aperture ground- 
plane. 


(*? - k l) |4i e ~ jkzi dl + elkzi dl ] - KK [/* dl + e jk * di ] = 

(*o - k 2 )T^ 0 e~ lkz ° dl - ktkyT^e-*'^ . (2.101) 


(2.43), (2.44), (2.87), (2.88), (2.100), and (2.101) form a linear system of six 
equations in six unknown spectra, given by the matrix equation: 


*11 

0 

-*11 

0 

0 

0 




v 

0 

-*11 

0 

*11 

0 

0 


I<j>l 


0 

*31 

0 

*33 

0 

*35 

0 




0 

0 

*31 

0 

*33 

0 

*35 




0 

*51 

*52 

*53 

*54 

*55 

*56 




0 

. *61 

*51 

*63 

*53 

*65 

*55 


** 


0 


( 2 . 102 ) 
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where 


*11 

= -jK » 

254 

1 

<N rH 

II 

k 2 x )e +jk ‘i dl , 


*31 

= —jk Zl e~* k ‘ l dl , 

Z 55 

= k x k y e 

-jk^di 


*33 

= jk Zl e + ’ k ‘> d ', 

2 5 6 

= -« 

— kl)e ~* k ‘° dl , 


*35 

= jK e ~ ikz ° dl » 

261 

= «- 

k 2 v )e-^ dl , 

(2.103) 

*51 

= -k x k y e~ jk ^ dl } 

263 

= (*?- 



*52 

= {k\-kl)e-^ 

265 

= -(K 

- k y )e~i kl ° dl . 


*53 

= -k x k v e+’ k " d \ 





matrix equation of (2.102) 

can 

be reduced 

into a two equation 

system in the 


two desired unknowns I^ x and . Thus the solution for the spectra /^ l and is 


given by 




V 1 Z 22 — ^2^12 

Z\\Z22 ~ Z\2 

(2.104) 

11 

V 2 Z 1 I — V\Z\2 
Z\\Z22 ~ Z \ 2 

(2.105) 

where 



Z u = 

(kj — k 2 )k z 0 cos(k I1 di) + j{ko ~ k x )k Zl sin(k Zl di), 

(2.106) 

Z\2 ~ 

-k x k y [k Zo cos(k Zl di) + jk Zl sin(fc Zl di)] , 

(2.107) 

Z 22 — 

(*1 - k 2 y )k zo cos(k z , di) + i(fco ~ k l) k *i ■“(*« i^)* 

(2.108) 

V, = 

k x k y (k zo + k z , ) 

(2.109) 

v 2 = 

[k Z0 (fc 2 - kj) + fc Zl (kg - kg)] 2k jf. 

(2.110) 


Note that the substitution of variables in (2.38) and (2.39) is needed in (2.106)- 
(2.110) for the purposes of computation of the aperture admittance. 


2.1.5 One Dielectric Layer - PEC Sheet Backed 

Figure 2.7 shows one dielectric layer, PEC sheet backed, on a rectangular aperture 
groundplane. The width of the single layer is d\. The incident spectra 7^, and 7^ 
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are determined by the enforcement of the boundary conditions at the two interfaces 
at z — 0 and d\. The boundary condition enforcement occurs at each set (fc x ,Jk y ) of 
spectral frequencies. The enforcement of boundary conditions results in a system of 
linear equations with the spectra as unknowns to be solved. The boundary conditions 
and resulting equations at z = 0 and are the same as for the free space backed case 
of 2.1.1. 



Figure 2.7: One dielectric layer, PEC sheet backed on a rectangular aperture ground- 
plane. 


The boundary condition at z = d x is that of vanishing tangential electric field, 
which produces the two equations: 


-j k Zl 7 0l e~* k ‘ l dl + jk Zl e jk *i dl = 0, 


( 2 . 111 ) 


jk Zl I^e lk ‘ ldl - j k 2l e jk *' d< = 0. 


( 2 . 112 ) 
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(2.43), (2.44), (2.111), and (2.112) form a linear system of four equations in four 
unknown spectra, given by the matrix equation: 


where 


*11 

0 

0 

“*11 

*31 

0 

. 0 

■231 

*11 

= ~jk. 

*31 

= ~jk. 

*33 

= iK 1 


-*n 

0 

*33 

0 


*1 > 


0 

■211 

0 

233 



■ 




Ifa 


0 


^1 


0 




0 


The solution of (2.113) for the desired spectra and 1 ^ is 

e +jk Zl di j 

^ 2k Zl sin(k Zl diY 

I<n = 0 . 


(2.113) 


(2.114) 


(2.115) 

(2.116) 


Note that the substitution of variables in (2.38) and (2.39) is needed in (2.115) for 
the purposes of computation of the aperture admittance. 
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2.2 Coaxial Aperture 


The admittance of a coaxial aperture is developed here. The analysis involves rep- 
resenting the fields external to the aperture in terms of the Fourier integral of their 
plane wave spectra. The aperture admittance is obtained by using a variational 
assumption for the aperture total field distribution and matching the boundary con- 
ditions for the various dielectric slab/resistive sheet configurations under study. 

The analysis starts with the representation of the electric vector potential F in 
terms of the cylindrical electric scalar potential ip: 

F = M- (2.117) 

It will become evident that since the fields in the coaxial aperture are axially sym- 
metric for the purposes of this investigation, only the <j> component of F is required. 
Referring to Figure 2.8, it can be seen that the potential is a functions of p and 
z, defined only for z > 0. The subscript 1 refers to Region 1, the dielectric slab 
in contact with the groundplane. The subscript 2 refers to Region 2, a second di- 
electric slab which may be on the first dielectric slab. The free space region above 
the dielectric slab(s) is denoted as Region 0 and all quantities associated with that 
region possess the subscript 0. The scalar potential is a solution to the scalar wave 
equation 

(V + ^ 0 , 1 , 2 )^ = 0, (2.118) 

where the wavenumbers &o,i,2 are given by 

k 2 0 = w 2 /ioeo, (2.119) 

for Region 0, and 

*1% = “>Voei,2 = *0 e r li2 , (2.120) 

for Regions 1 and 2. In general the relative permittivities e ri are complex, repre- 
senting possible losses in the slabs. 
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REGION 1 (€ 1 ,fl 0 ) 

POSSIBLE RESISTIVE 
OR PEC SHEET 

REGION 0 (e 0 ,/U q) 



REGION 1 { e 1 , fJ, Q ) 
REGION 2 {£ 2 ,/J. c ) 
POSSIBLE PEC SHEET 


REGION 0 


Figure 2.8: Dielectric slab/resistive sheet configurations on coaxial aperture ground- 
plane. 
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The scalar potential can be expressed in terms of its plane-wave spectra. In 
Regions 1 and 2, the total spectra is comprised of the incident and reflected spectra, 
respectively denoted by / and R. In Region 0, the total spectra is comprised of one 
transmitted spectrum, denoted by T. The plane wave spectral representation chosen 
for the coaxial geometry involves the Fourier- Bessel transform with Bessel function 


of order 1. This representation has been chosen because J\(f3p) - 

0 as p 0, 

which is consistent with the fact that the fields vanish in the inner conductor of 
the coaxial aperture. This spectral representation is also consistent with the coaxial 
aperture work of Berrie [11]. Thus the plane wave spectral representation of the 

scalar potential in Region 1 is given by 


foe r 

Mp,*) = J o [hm^ z + R*APV kl ' z ] MPpWP-, 

(2.121) 

in Region 2 by 


foe r 

Mp, A = j o [ + RMe^] M0p)(3d/3; 

(2.122) 

and in Region 0 by 


MP,A = jH [T*(0)e-*- o'] MPp)pd(3. 

(2.123) 

The propagation constants are given by 


^*o,i,2 = ±\/*0,1.2 — 

(2.124) 

where the sign is chosen so that the radiation condition is satisfied, 
for each propagation constant is chosen so that 

Thus the sign 


(2.125) 

and 


S(ko,i,*) ^ °- 

(2.126) 

From vector potential theory, the electric and magnetic fields are 
electric vector potential by 

related to the 

E = -V x F = ^-p 

oz 

(2.127) 
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and 


H = t— — [fc 2 F + V(V • F)1 = (2.128) 

1 J 

Therefore the non-zero cylindrical components of the electric and magnetic fields in 
Region 1 are given by 

E„(p,z) = jH + jfc, «*,€»-*] J,(0p)0d0, (2.129) 

and 

/f<tf>i(^J z ) = J ju;ei Ji(/3p)/3d/3; (2.130) 

in Region 2 by 

£„(*>. z) = f + }k,,Rp,e*i‘] M0p)0d0, (2.131) 

and 

H^faz) = J™ [-jvc2l*e- ik *‘ - M0p)0d0; (2.132) 

and in Region 0 by 

E„(p, z) = jH M0p)0d0, (2-133) 

and 

H^(p, z) = jH [-ju,€„T* I e-'‘- < ] M0p)0d0. (2.134) 

Note that in the above the 0 functional dependence is dropped from the spectra and 
is understood. 

As in the rectangular aperture case, the solution to the fields can be obtained 
by satisfying the boundary conditions at the groundplane/dielectric slab interface 
and other interfaces present in the dielectric slab/resistive sheet configurations, but 
only if the total fields in the aperture are known. However, the total fields in the 
aperture are not known, since they consist of an infinite sum of waveguide modes with 
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unknown coefficients. In order to find the coefficients, a Moment Method solution 
can be set up in a manner similar to that which has been performed for the coaxial 
aperture case. For the purposes of this investigation a variational solution will be 
used [3] and its results compared to the Moment Method solution of Berrie [11]. The 
variational solution assumption for the coaxial aperture total electric field is of the 
form of the dominant TEM mode. The importance of higher-order coaxial modes is 
investigated in Subsection 3.3.2. 

Thus the total electric field in the plane z = 0 is assumed to be 


E(p>0) — e o — < 


where the factor 


in the aperture 


\/2w ln(b/a)p ’ 

0 , outside the aperture 

normalizes the electric field so that 

1. 


(2.135) 


Y/2?rIn(6/a) 

i»27T [OQ 

/ / E(p, 0) - E(p, 0 )pdpd(j> = 1. (2.136) 

Jo Jo 

Referring back to Figure 2.8, it can be seen that a and b are the inner and outer 
radii of the coaxial waveguide, respectively. 

The coaxial waveguide aperture admittance, Y ai is given by 

/o Vo“ E(p>0) X H(p,0)M># 


Y„ = 


fOC 

= 2./ E„{p,0)n tl ( p< 0)pdp. (2.137) 

Jo 


[/o' /o“ E(p, 0) - 

Evaluation of (2.137) is accomplished through the use of Parseval’s Theorem in two 
dimensions in a cylindrical coordinate system independent of <f>: 


fOC #oo 

/ EpiH<tnpdp= S Pl 7ifal3d/3, 

Jo Jo 


(2.138) 


where and Ufa are the Fourier-Bessel transforms of E Pl and Hfa , at z = 0. These 
transforms can be found from (2.129) and (2.130) and are given by 


Zpi(ft) ~ [ dfa + 


(2.139) 


and 


(P) = -jwe i [I fa + J2*] 


(2.140) 
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Thus the aperture admittance can be written as 


Y a = 27T f°° {jk Zl [-J* + R* J} {—ju>ei [J* + i4,]} 0dj3. (2.141) 

JO 

The unknown quantities in (2.141) are the incident and reflected spectra of Region 
1. The boundary condition of continuous tangential electric field at z = 0 for each 
spectral value (3 gives the relationship between the Region 1 incident and reflected 
spectra as 


j^zi [ ] — 


MPa)- JoW 

yjfyrki (6/a)/? 


(2.142) 


where the Fourier-Bessel transform of the aperture electric field is denoted by / for 
the sake of brevity. Substitution of(2.142) into (2.141) allows the aperture admit- 
tance to be written as 


ln(6/a) Jo 


°° (-ja/e, [!<,, + J^]) (JoOJa) - MPb)? 


d0. 


(2.143) 


ln(6/a) Jo f /3 

In order to use the aperture admittance in the calculation of the reflection coef- 
ficient, it must first be normalized with respect to the dominant mode waveguide 
admittance: 


Y Y 
V = — = — 

" Y g Yo^Tc’ 


(2.144) 


where lo is the admittance in the waveguide and e rc is the relative permittivity of 
the waveguide material. Thus, the reflection coefficient of the coaxial waveguide 
aperture is given by 


r - rrf • < 2 - 145 > 

A few items about the integration in (2.143) must be stated. The (3 integration 
is performed along the path in the complex /3 plane shown in Figure 2.9. The 
integration path is deformed from the real axis so as to avoid surface wave poles. 
Surface wave poles can occur on the real axis when there is at least one lossless 
dielectric layer in any of the five configurations of Figure 2.1. These poles occur 
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in the interval kg < fl < e r ko, where e r is the greater of the two real relative 
permittivities, or is the only real relative permittivity. Loss in the dielectric layer(s) 
causes the layer(s) surface wave poles to move below the real axis. The integration 
path starts at 0 = 0, goes to ko + jOAko, then to |e rfl |*o + jOAko , 2|e rff |Jbo, and 
then oo + j 0. This integration path was determined experimentally and found to 
work well. The integration path crosses no poles during its deformation from the 
real axis, so no residues need to be computed and included in the integral. 


Im{p) 



Figure 2.9: Integration path in the complex (3 plane, with surface wave poles and 
parallel plate waveguide pole (PEC sheet backed cases only). 

If there is a PEC sheet present, an additional complication occurs when the 
frequency is high enough to excite one or more parallel plate waveguide modes. A 
mode manifests itself as a pole in the complex /3 plane by appearing at /? = 0 at 
the cutoff frequency, and then moves along the positive real axis as the frequency 
increases. No provision has been included in the codes for this contingency, since 
the codes are to be used only at frequencies where there are no propagating parallel 
plate waveguide modes. 

However, the solution for the reflection coefficient is by no means complete. 
The aperture admittance is still dependent on the unknown spectra and R t p 1 , 
which are dependent on the boundary conditions given by the particular dielectric 
slab/resistive or PEC sheet configuration in contact with the aperture ground plane. 
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The following subsections give the derivation of / ^ and R^ >1 for the five configurations 
under investigation. 

2.2.1 Two Dielectric Layers - Free Space Backed 

Figure 2.10 shows two dielectric layers, free space backed, on a coaxial aperture 
groundplane. The width of the first layer is d x , and the width of the second layer is 
d 2 — di. The spectra 7^ and are determined by the enforcement of the boundary 
conditions at the three interfaces at z = 0, d\, and d 2 . The boundary condition en- 
forcement occurs at every spectral value /?. The enforcement of boundary conditions 
results in a system of linear equations with the spectra as unknowns to be solved. 
The boundary conditions with their resulting equations are listed here: 

The boundary condition of continuous tangential electric field across the interface 
at z = 0 gives the equation: 

=/• (2.146) 
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The boundary condition of continuous tangential electric field across the interface 
at z = d\ gives the equation: 

-k z ,R^e jh 'i d ' = - k^R^e*”* . (2.147) 

The boundary condition of continuous tangential magnetic field across the interface 
at z = di gives the equation: 

*1 /*«-*■** + e Tl R^e jk 'i dl = e r2 I^ e - jk ‘* dl + e ri R^e jk ’* dl . (2.148) 

The boundary condition of continuous tangential electric field across the interface 
at z = d 2 gives the equation: 

k Z7 I xh e~ ik " d ' - k Z7 R^e jk ‘> d * = k Zo T^ 0 e~ jk ‘o d \ (2.149) 

The boundary condition of continuous tangential magnetic field across the interface 
at z = d 2 gives the equation: 


e r7 I^e- jk ‘* d * + e r7 R« J7 e jk '* d > = T^e~ jk ‘° d \ 


(2.150) 


These five linear equations in the five unknown 
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spectra form this matrix equation: 


(2.151) 


where 
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= 
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e r7 e~ jk ^ d \ 

z 54 
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£ p+i**2^2 

t r2 C , 

z 55 

— 

—e~ 3k, ° dl . 


(2.152) 
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In the code (2.151) is solved numerically for the spectra I ^ and R ^ 1 . 

2.2.2 Two Dielectric Layers - PEC Sheet Backed 

Figure 2.11 shows two dielectric layers, PEC sheet backed, on a coaxial aperture 
groundplane. The width of the first layer is di, and the width of the second layer is 
d 2 -d\. The PEC sheet is located at z = d 2 . The spectra 1^ and R^, l are determined 
by the enforcement of the boundary conditions at the three interfaces at z — 0, dj, 
and d 2 . The boundary condition enforcement occurs at each spectral value /3. The 
enforcement of boundary conditions results in a system of linear equations with the 
spectra as unknowns to be solved. The boundary conditions and resulting equations 
at z = 0 and z — d\ are the same as for the free space backed case of 2.2.1. 
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(2.146)-(2.148), and (2.153) form a linear system of four equations in four un- 
known spectra, given by the matrix equation: 

2u Z\2 0 0 

z 2l z 22 z 23 z 24 

Z 31 z 32 z 33 z 34 

0 0 243 244 

where 

z u = -jk z ,, 
z n = jk Zl , 
z 2i = k Zl e- jk *i dl , 
z 22 = -k Zl e +jk ‘i dx 
z 23 = -k Z7 e~ jk ^ dl 
2 24 = k Z2 e +jk ** dl , 

In the code (2.154) is solved numerically for the spectra 1^ and R^ l . 

2.2.3 One Dielectric Layer - Resistive Sheet Backed 

Figure 2.12 shows one dielectric layer, resistive sheet backed, on a coaxial aperture 
groundplane. The width of the single layer is d\ . The resistive sheet is located at 2 
= di . The spectra and are determined by the enforcement of the boundary 
conditions at the two interfaces at 2 = 0 and d\ . The boundary condition enforcement 
occurs at each spectral value /3. The enforcement of boundary conditions results in a 
system of linear equations with the spectra as unknowns to be solved. The boundary 
conditions and resulting equations at 2 = 0 and are the same as for the free space 
backed case of 2.2.1. 

The resistive boundary conditions that the electric and magnetic fields obey at 
2 = d\ are 
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— 1 


R^> 1 


0 


ifa 


0 


r ^ 2 


0 


Z 31 = 

f p-jkz ! di 

c ri c ? 

z 32 = 

e ri e +jk ‘i d 1 , 

z 33 = 

— e ri e~* kl ? dl 

z 34 = 

—e r7 e + i k, ? dl 

z 43 = 

k Z7 e ~ jk ** d * , 

z 44 — 

-k Z2 e +]k ^ d2 


(2.154) 


(2.155) 


z x (E 0 - Ei) = 0 


(2.156) 
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REGION 1 {€ 1 ,fjL 0 ) 
RESISTIVE SHEET 

REGION 0 ( € q , /^q) 

Figure 2.12: One dielectric layer, resistive sheet backed on a coaxial aperture ground- 
plane. 

z x z x E = —Ri x (H 0 — Hi) (2.157) 

(2.156) gives rise to this equation: 

KhS~ ik " dx ~ KRo^ k " dl = (2.158) 

(2.157) gives rise to this equation: 

vdRI* e- jk ‘' dl + e +jk " dl = [u>eo R + **] T^e~ ik ^ dl . (2.159) 

(2.146), (2.158), and (2.159) form a linear system of three equations in three 
unknown spectra, given by the matrix equation: 


-jk z , jk Zl 0 
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f 

k Xi e-’ k ‘i dl —k Zl —k Zo e~* kl ° dl 



= 

0 

weiRe- ik * dl u>eiRe + * k * 1 dl - [ue 0 R + k Zo ] e~’ kt * dl 


i 

i 


0 


In the code (2.160) is solved numerically for the spectra 7^, and . 
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2.2.4 One Dielectric Layer - Free Space Backed 

Figure 2.13 shows one dielectric layer, free space backed, on a coaxial aperture 
groundplane. The width of the single layer is d ] . The matrix equation for the three 
spectra R ^, l , and T^, 0 is readily obtained by letting the resistivity R in (2.160) 
tend to oo, giving: 


jk Zl jk Zl 0 


1 

i 


» 

k Zl e~ jkx ' dl - k Zl e jk *i dl -k Zo e~ jkx o dl 



ii 

0 

ti e ~* kz i dl e +J,Ar *i dl — e 0 e~^ kz o dl 


i 

i 


o 

1 


In the code (2.161) is solved numerically for the spectra I ^ and R^ x . 


(2.161) 



Figure 2.13: One dielectric layer, free space backed on a coaxial aperture ground- 
plane. 
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2.2.5 One Dielectric Layer - PEC Sheet Backed 


Figure 2.14 shows one dielectric layer, PEC sheet backed, on a coaxial aperture 
groundplane. The width of the single layer is d \ . Since only the sum in 

(2.143) is necessary for the admittance calculation, and is in a compact form readily 
obtained by letting R go to 0 in (2.160), it is presented here: 


I+ + R* tan(fc 2l di) 


(2.162) 


This spectral sum is the expression employed in the code for the single PEC sheet 
backed dielectric layer case. 




d 
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REGION 1 {£ 1 ,fJ, 0 ) 
PEC SHEET 


Figure 2.14: One dielectric layer, PEC sheet backed on a coaxial aperture ground- 
plane. 
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Chapter 3 

Parameter Study 


A parameter study employing the numerical model of the previous chapter has been 
conducted. The parameter study focuses on three topics. The first is the behavior of 
the reflection coefficient for a single free space or PEC sheet backed dielectric layer 
when permittivity is a parameter fixed over frequency. The second is the validity of 
the approximation of a thin dielectric layer by a thin resistive sheet. The third is a 
study of the errors introduced by air gaps and the neglect of higher order modes. 

3.1 Single Layer 

In this section the single layer geometries of Figures 2.6, 2.7, 2.13, and 2.14 are 
studied. In all cases the thickness of the dielectric sheet is di = .2". The permittivity 
parameter e r is set at either 2, 4, 6, or 8. The rectangular aperture employed is a free 
space filled X band waveguide aperture with a = A" and b = .9" at / = 8-12 GHz. 
The coaxial aperture employed is a Teflon (e r = 2.1) filled 50ft coaxial waveguide 
aperture with a = .5", b = 1.153" at / = .1-2 GHz. 

Figure 3.1 shows the reflection coefficient for the rectangular aperture, free space 
backed case. Figure 3.2 shows the reflection coefficient for the rectangular aperture, 
PEC sheet backed case. These figures illustrate that the reflection coefficient show a 
strong dependency upon the permittivity. This is also a very desirable feature since 
the extraction process will not be as sensitive to measurement noise. The sample 
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thickness determination is based upon two criteria. The first criterion is to choose a 
thickness where no nulls are present. The second criterion is to limit the thickness 
to less than a quarter wavelength in the sample to avoid the generation of surface 
waves. The presence of surface waves is detrimental for a finite groundplane/sample 
size. This is due to the surface waves being reflected back to the aperture which 
changes the aperture impedance from the value for an infinite coated groundplane. 

Figure 3.3 shows the reflection coefficient for the coaxial aperture, free space 
backed case. Figure 3.4 shows the reflection coefficient for the coaxial aperture, 
PEC sheet backed case. The dependency of the permittivity upon the reflection 
coefficient is not very strong for the free space back case, at least for this frequency 
range and coaxial aperture size. More desirable results are obtain when the sample 
is PEC backed. 
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Figure 3.3: Teflon filled 500 coaxial aperture, .2" slab, free space backed, 
2-solid, 4-dotted, 6-dashed, 8-dotdashed 
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Phase (degrees) 
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Figure 3.4: Teflon filled 5011 coaxial aperture, .2" slab, PEC sheet backed, e r = 
2-solid, 4-dotted, 6-dashed, 8-dotdashed 
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3.2 Dual Layer 


In this section the resistive sheet backed single layer geometry of Figure 2.5 is com- 
pared with the free space backed dual layer geometry of Figure 2.3. The aperture 
employed is a free space filled X band rectangular waveguide aperture with a = .4" 
and b = .9". The frequency range is / = 8-12 GHz. The dielectric layer in contact 
with the groundplane has a thickness of .2" and a relative permittivity of e r = 4. 
The frequency range is / = 8-12 GHz. The resistive sheet in contact with the single 
dielectric layer has a resistivity of R = 10 3 fi/D, 10 4 fI/D, or 10 5 fi/n. The second 
dielectric layer in the dual layer structure is made of conducting material only;, i.e., 



(3.1) 


where a is the conductivity of the material in Siemens (fi _1 ). The thickness of the 
second layer is t — .002". The conductivity of the second layer is chosen to satisfy 
the resistive sheet approximation to the thin material slab, utilizing the resistivity 
of the resistive sheet in the single layer configuration: 



(3.2) 


which holds as long as t satisfies the condition [13] 

l*</2| = « 1- 


(3.3) 


Thus the three values of a corresponding in order to the respective R values above 
are 19.7 1.97 fl -1 , and .197 fl -1 . The quantity \kt/2\ for each of the three 

values of a is plotted versus frequency in Figure 3.5. 
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Figure 3.5: \kt/2\ versus frequency - c = 19.712 1 - solid, <r = 1.9712 1 - dotted, a 
— .19712“ 1 - dashed. 

Figure 3.6 shows the reflection coefficient for the single layer, resistive sheet 
backed cases. Figure 3.7 shows the reflection coefficient for the corresponding dual 
layer cases. Note that for the case where |fct/2| is the smallest, (<r = .19712 -1 ) the 
thin sheet approximation holds. The thin sheet approximation also performs well for 
the middle <r value, but not well for the largest a value, where \kt/2\ is the largest, 
on the order of 0.3. This demonstrates the necessity of satisfying the condition that 
\kt/2\ must be much less than 1 in order for the thin sheet approximation to be 
applicable. Note for a thicker sample, say .02" thick, the range of R values can vary 
from R = 10 2 to 10 4 12/D to have the same performance as the shown cases. 

The dependency of the reflection coefficient upon the resistivity of a resistive 
sheet is demonstrated in Figure 3.8. The resistive sheet varies in resistivity from 
012/ □ to 300012/ □, and is in contact with a .2" thick dielectric sheet with e r = 
4. The dielectric sheet rests on an X-band rectangular waveguide aperture with a 
dominant mode excitation of 10 GHz. 
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Phase (degrees) 


Frequency (GHz) 



Figure 3.6: Free space filled rectangular X band aperture, Resistive sheet backed, 
•2" slab with e r = 4, R = 10 3 ft/n - solid, R = 10 4 fi/D - dotted, R = 10 s fi/D - 
dashed. 





Phase (degrees) Magnitude (dB) 




Figure 3.7: Ftee space filled rectangular X band aperture, .002" sheet backed, .2" 
slab with e r = 4, cr = 19. 7fi~ 1 - solid, <r = 1.97fi~* - dotted, <r = .197fi 1 - dashed. 
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Figure 3.8: Reflection coefficient dependency on resistivity at 10 GHz for a resis- 
tive sheet backed .2" thick dielectric sheet with e r = 4 on an X-band rectangular 
waveguide aperture. 
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3.3 Errors 


In this section the errors introduced by air gaps and the neglect of higher order modes 
are studied. The studies of these two error sources are of practical value. The air 
gap study will examine the sensitivity of the reflection coefficient measurement to 
uniform air gaps between the groundplane and the dielectric layer that should be 
flush with the groundplane. The higher order mode study will examine the limits in 
frequency where the lowest order mode approximation is valid. 

3.3.1 Air Gaps 

The air gap study involves the single layer free space backed rectangular aperture 
geometry of Figure 2.6. The rectangular aperture employed is a free space filled 
rectangular X band waveguide aperture with a = A" and b = .9" at / = 8-12 GHz. 
The coaxial aperture employed is a Teflon filled 50fl coaxial waveguide aperture with 
a = .5", 6 = 1.153" at / = .1-2 GHz. The dielectric sheet has a thickness of d\ = 
.2" and relative permittivity of either e r = 4 or 8. 

The reflection coefficient for the rectangular aperture case is shown in Figures 3.9 
and 3.10 for e T = 4 and 8 respectively. The reflection coefficient with no air gap is 
compared against the reflection coefficients for gaps .05% and .10% the thickness of 
the dielectric layer. Note that in Figure 3.10 magnitude information near the null is 
not shown since the null does not contain any useful information. The permittivities 
extracted from the reflection coefficient for the gap cases assuming there is no gap 
are shown in Figures 3.11 and 3.12 for e r = 4 and 8 respectively. It can be seen that 
as the gap width is increased the extraction becomes more inaccurate. Also, the gap 
has a greater effect for higher permittivities. Therefore, a higher permittivity sample 
requires more care in the prevention of air gaps to ensure an accurate permittivity 
extraction. 

The reflection coefficient for the coaxial aperture case is shown for e r = 4 and 8 
in Figures 3.13 and 3.14, respectively. The reflection coefficient with no air gap is 
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compared against the reflection coefficients for gaps .05% and .10% the thickness of 
the dielectric layer. The permittivities extracted from the reflection coefficient for 
the gap cases assuming there is no gap are shown in Figures 3.15 and 3.16 for e r = 4 
and 8 respectively. Note that the errors due to the air gap are not manifested in the 
real part of the permittivity, but in the imaginary part. As the gap width increases, 
the error increases. The error is greater for the greater permittivity. Note that great 
care must be taken to avoid air gaps for the coaxial measurement configuration for 
even low permittivity samples, since even small gaps produce a significant extraction 
error. 
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Frequency (GHz) 

Figure 3.9: Air gap study, X band, .2", e r = 4. No gap-solid, .05% gap-dotted, .10% 
gap-dashed. 







Frequency (GHz) 


Figure 3.10: Air gap study, X band, .2", e r = 8. No gap-solid, .05% gap-dotted, 
.10% gap-dashed. 
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Frequency (GHz) 


Figure 3.12: Permittivity extraction, X band, .2", e r = 8. No gap-solid, .05% gap 
dotted, .10% gap-dashed. 
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Figure 3.15: Permittivity extraction, coaxial, .2", e r = 4. No gap-solid, .05% gap 
dotted, .10% gap-dashed. 
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Figure 3.16: Permittivity extraction, coaxial, .2", e r = 8. No gap-solid, .05% gap- 
dotted, .10% gap-dashed. 
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3.3.2 Higher Order Modes 

The existence of higher order modes can occur for both the rectangular and coaxial 
apertures. Their effect was observed during the calibration of the measured data for 
the rectangular guide. Generally, three measurements are required to calibrate data. 
The three measurements originally taken were a short, open and matched load at the 
aperture. The resulting calibrated data typically had an oscillation which was not 
explainable. A potential source of error existed in the calibration due to the “exact” 
reference used for the open measurement. The exact reference value assumed only 
a TEio in the waveguide at the aperture. The approximation is very reasonable 
and sufficiently accurate for many applications. However, for inverse problems (such 
as material parameter determination), this may not be sufficient as shown in the 
following figures. 

The calculated (single mode only) and measured reflection coefficient in the X- 
band is shown in Figure 3.17. Note that the measured response is offset and lower 
than the calculated response by approximately 1 dB. An additional calculation us- 
ing a hybrid FEM/BEM code for an aperture in an infinite groundplane at 10 GHz 
yielded a reflection coefficient of -12.8 dB which is in agreement with the measure- 
ment. The calibration of the measured reflection coefficient for the open aperture 
used three calibration standards of a short, an offset short and a matched load. Fig- 
ure 3.18 illustrates the effect of calibrating the measured reflection coefficient for a 
dielectric (lossy foam, thickness = .595 in.) covered aperture with two sets of cali- 
bration standards. One standard set uses a short, an offset short and a matched load 
where the other standard replaces the offset short with an open aperture (free space) 
measurement. It is clearly seen that the calibration with the open aperture mea- 
surement has a non-physical oscillation. This oscillation can be directly attributed 
to using only one expansion mode in the waveguide for the exact reference value for 
the open aperture measurement. 
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Figure 3.17: Reflection coefficient for a X-band aperture. Measured - solid, calcula- 
tion (single mode expansion) - dashed, hybrid FEM/BEM calculation - star. 
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The necessity of higher order modes in the analytical solution is very apparent. 
Most references pertaining to this technique ignore the effect of higher order modes. 
Of two references that discuss the effect of higher order modes for rectangular aper- 
tures, only one has results that demonstrate their importance [8], especially when 
the aperture has a material covering. The other reference does not consider them 
to be very important [12]. It should be noted that the formulation and computer 
software for the coaxial aperture does include the effect of the higher order modes. 
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Chapter 4 
Measurements 


The measurement procedure and some results are presented in this chapter. In the 
first section, a brief discussion of the measurement configuration and calibration is 
given. In the second section, the Newton- Raphson extraction procedure is described. 
In the third section, some measured results are presented along with a comparison 
against the numerical model. 

4.1 Measurement Configuration and Calibration 

The measurement configuration consists of a Hewlett-Packard 8510B network ana- 
lyzer with a groundplane fixture having a rectangular aperture in the center of the 
groundplane. The groundplane measures 9" square to simulate an infinite ground- 
plane in the 8 to 12 GHz band (X-band). These dimensions allow the separation 
between the aperture to the groundplane edge to be several wavelengths. Addition- 
ally, the perimeters of the aperture and the groundplane are rotated 45 degrees to 
minimize any direct surface wave/edge diffraction back into the aperture. 

The reflection coefficient data must be calibrated before de-embedding of the 
permittivity can be performed. There are three types of linear errors in reflection 
coefficient measurements: directivity, frequency response, and source match, rep- 
resented respectively by the dimensionless quantities E D , Er and E s . All three 
error terms are complex and frequency dependent, and ideally Ed and Es are 0 and 
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Er is 1. The effect of these three linear errors on the actual reflection coefficient 
Sn a , which produces the measured reflection coefficient Sn m , can be modeled as the 
two-port error adapter of Figure 4.1. 



Figure 4.1: One port error model 


The values of Ep, Er, and Es must be determined for each frequency so that 
the data can be calibrated. Sn a is related to Sn m through this equation: 

Ep + ' Efi5ll<1 = S n (4.1) 

Three such equations as in (4.1) are required for solution of the three unknowns Ep, 
Er, and Es- These are formed by taking three reflection coefficient measurements 
for which Sn a is known. The three reflection coefficient measurements employed 
in the calibration code developed are the short circuit, offset short circuit, and 
matched load. The short circuit is formed by placing a metal plate against the 
waveguide aperture. The metal plate is typically aluminum and considered PEC so 
that Sn 0 = —1. The offset short circuit is formed by extending the waveguide 1.2 
cm through the use of a material sample holder and shorting the open end. In this 
case Su a = —e~i k2d where d = .012 m. The matched load is formed by attaching to 
the aperture a lossy waveguide termination, with an assumed reflection coefficient 
S Ua = 0. With the error terms determined, they are used in the solution of (4.1) 


for S Ua - 

S = Slim — E D 

Es(Sn m -E D ) + E R ’ 

This completes the calibration of the reflection coefficient data. 
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4.2 Newton-Raphson Extraction Procedure 


The permittivity of a material sample is extracted from its calibrated reflection 
coefficient data via the Newton-Raphson iterative procedure for finding roots of 
transcendental equations [15]. The general form of the iterative Newton-Raphson 
equations are 


/(*) = o, 


(4.3) 


and 


^n+1 — 


fM 


(4.4) 


where / is a function of x. It is desired to find an x that satisfies (4.3), which is 
called a root. The Newton-Raphson procedure is initiated with a guess of the root, 
called Xo. xo is inserted into (4.4), which produces x\. The process by which x n 
is obtained is called the n th iteration. If the function is properly behaved, the x n ’s 
would converge to a zero of (4.3). A convergence criterion on x n is employed to 
terminate the Newton-Raphson method when the desired level of accuracy for the 
zero is achieved. 

The permittivity of the material sample is defined to be the zero of this functional: 


f{e h ) = T(e = e h )-T(e=e g ) (4.5) 

where T is the reflection coefficient of the aperture at a certain frequency, which 
is covered by a dielectric slab of permittivity e and the thickness of the material 
sample slab. e h is the permittivity of a homogeneous slab of the same thickness as 
the sample slab of unknown permittivity e„. 

The use of T can be incorporated in one of three ways in (4.5). The first is just 
to use the magnitude of T which is appropriate for a lossless dielectric. The phase 
of r can also be used when the dielectric is lossless. However, when the dielectric is 
lossy, then the complex value of T is required. 
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A c/, is sought via the Newton- Raphson procedure so that 

/(efc) = 0 


(4.6) 


is true. In order to apply the Newton- Raphson procedure to (4.6), the derivative of 
f(th), must be computed. It is impractical to compute f'(ch) exactly, however, 

a simple approximation to /^(e/,) is given by the following central difference formula: 

f( e h) « /(e, ' + Ae) 2 ^ /(£ ''~ — ' W 

where Ae is a differential parameter such that |Ae| *C 1. 

The convergence criterion of the Newton- Raphson procedure is given by 

i/(«)i < e < 4 - 8 ) 

where e is the convergence criterion constant, a small real number such that |e| 1* 

At the end of the first iteration when the convergence criterion is met, the Newton- 
Raphson procedure is terminated and e h is called the permittivity of the sample 
material. 


4.3 Examples 

The per mi ttivity of three samples were extracted from calibrated reflection mea- 
surements in the X-band. Prior to de-embedding, the calibrated data was smoothed 
using a moving average to minimize the effect of noise. The samples were acrylic, 
epoxy resin, and a carbon loaded foam. The dimensions of the samples are shown in 
Table 4.1. It should be noted that the accuracy of the material value de-embedding 
is limited due to the use of a single waveguide mode expansion. 

The first example is that of the permittivity extraction of a .111 in. thick acrylic 
plexiglass sheet. A nominal relative permittivity value for plexiglass at 10 GHz is 
e r = 2.59 — j'0.017 [14], and varies insignificantly from this value throughout the 
X band. The raw and calibrated reflection coefficient data is presented in Figure 
4.2, along with a comparison against the reference solution for a dielectric sheet of 
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the same thickness and permittivity. Note that a significant difference occurs in the 
magnitude between the calibrated and calculated (using the nominal permittivity 
value) reflection coefficient. One source of error can be due to an air gap between 
the sample and the groundplane. However, the influence of higher order modes can 
also account for this and has been shown to be important. 

The de-embedded relative permittivity is shown in Figure 4.3 where the complex 
reflection coefficient was used to obtain both the real and imaginary part of the 
permittivity. The average de-embedded relative permittivity is approximately 2.3 - 
jO.l. The discrepancy between the de-embedded and nominal relative permittivities 
is most likely due to the neglect of higher order modes. 

The second example is that of the permittivity extraction from a .256 in. thick 
epoxy resin sample. Figure 4.4 shows the epoxy’s calibrated reflection coefficient 
data compared against reference curves of the reflection coefficient for a dielectric 
slab .256 in. with permittivity e r = 8.5. Note that the calibrated reflection coefficient 
is corrupted in the region 8-9 GHz, where there is a minimum in the magnitude of 
the reflection coefficient. The extracted permittivity of the epoxy is shown in Figure 
4.5. The permittivity for frequencies below 8.5 GHz were not calculated due to the 
difficulties introduced by the corruption. Note that the epoxy is found to be basically 
lossless with a relative permittivity ranging from 9.5 at 8.5 GHz down to 8.0 at 12 
GHz. Note that the imaginary component of the dielectric constant is positive for 
some frequencies. This behavior can be attributed to error introduced when higher 
order modes are ignored in the aperture. The effect of these the higher order modes 
become more apparent for electrically thick material coatings [8]. 


71 



Phase (degrees) Magnitude (dB) 




Figure 4.2: Reflection coefficient data for an X band rectangular aperture reflection 
measurement of a plexiglass sheet .111 in. thick. (Raw - solid, calibrated - dotted, 
reference - dashed.) 
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Phase (degrees) Magnitude (dB) 


0 



Figure 4.4: Reflection coefficient data for an X band rectangular aperture reflection 
measurement of an epoxy resin sheet .256 in. thick. (Raw - solid, reference - dotted.) 
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Table 4.1: Material Sample Dimensions 


Material 

Thickness (inches) 

Sides (inches) 

Acrylic 

0.111 


Epoxy 

0.256 

4.9 x 5.0 

Lossy Foam 1 

0.595 

12 x 14 

Lossy Foam 2 

0.886 

12 x 14 
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The third example is that of the permittivity extraction from two samples of a 
commercially available lossy foam (anechoic chamber foam), one .595 in. thick, the 
other .886 in. thick. The relative permittivity of the lossy foam obtained via slotted 
line is approximately e r = 1.2 - jO.38. This value is extrapolated from measurements 
between .1 to 6 GHz. Figure 4.6 shows the calibrated reflection coefficient data 
for the lossy foam samples and Figure 4.7 shows the permittivities extracted from 
the complex reflection coefficient of the lossy foam samples. Similar de-embedded 
permittivity values were obtained for both sample thicknesses. The average de- 
embedded relative permittivity is approximately 1.05 - j0.175. The discrepancy 
between the de-embedded and slotted line method relative permittivities is most 
likely due to the neglect of higher order modes. 
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Phase (degrees) Magnitude (dB) 



Frequency (GHz) 

Figure 4.6: Reflection coefficient data for an X band rectangular aperture reflection 
measurement of lossy foam sheets. (.595 in. thick - solid, .886 in. thick • dotted.) 
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lm{Relative Permittivity} Re{Relative Permittivity} 




Figure 4.7: Permittivity extraction from two lossy foam sheets via X-Band rectan- 
gular aperture. (.595 in. thick - solid, .886 in. thick - dotted.) 
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Chapter 5 


Conclusions 


A variety of analytical reflection coefficient solutions were presented for material 
loaded rectangular apertures in an infinite groundplane. These solutions provide an 
analytical reference in a Newton-Raphson based technique to extract the material 
value from the material. The analytical solutions, for simplicity, involved only the 
dominate TE 10 mode in the waveguide. The restriction limits the accuracy of the 
resulting extracted values. The inclusion of the higher order modes decreases the 
attractiveness of this general approach do to a significant computational requirement. 

However, some useful material value results Me still possible considering the 
simplicity of the Newton-Raphson technique and its non-destructive nature to planar 
material samples. 
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